Soliton absorption spectroscopy 
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We analyze optical soliton propagation in the presence of weak absorption lines with much nar- 
rower linewidths as compared to the soliton spectrum width using the novel perturbation analysis 
technique based on an integral representation in the spectral domain. The stable soliton acquires 
spectral modulation that follows the associated index of refraction of the absorber. The model can 
be applied to ordinary soliton propagation and to an absorber inside a passively modelocked laser. 
In the latter case, a comparison with water vapor absorption in a femtosecond CriZnSe laser yields 
a very good agreement with experiment. Compared to the conventional absorption measurement in 
a cell of the same length, the signal is increased by an order of magnitude. The obtained analytical 
expressions allow further improving of the sensitivity and spectroscopic accuracy making the soliton 
absorption spectroscopy a promising novel measurement technique. 



1. INTRODUCTION 

Light sources based on femtosecond pulse oscilla- 
tors have now become widely used tools for ultrashort 
studies, optical metrology, and spectroscopy. Such 
sources combine broad smooth spectra with diffraction- 
limited brightness, which is especially important for high- 
sensitivity spectroscopic applications. Advances in near- 
and mid-infrared femtosecond oscillators made possible 
operation in the wavelength ranges of strong molecu- 
lar absorption, allowing direct measurement of important 
molecular gases with high resolution and good signal-to- 
noise ratio [if. At the same time, it was observed that 
such oscillators behave quite differently, when the absorb- 
ing gas fills the laser cavity or introduced after the output 
mirror 0, Q . The issue has become especially important 
with introduction of the mid-IR femtosecond oscillators 
such as CriZnSe Q, which operate in the 2-3 ^m wave- 
length region with strong atmospheric absorption. 

As an example, Figure [T] presents a typical spectrum 
of a Cr:ZnSe femtosecond oscillator, operating at normal 
atmospheric conditions. It is clearly seen, that the pulse 
spectrum acquires strong modulation features which re- 
semble the dispersion signatures of the atmospheric lines. 
Being undesirable for some applications, such spectral 
modulation might at the same time open up interest- 
ing opportunity of intracavity absorption spectroscopy. 
Compared with the traditional intracavity laser absorp- 
tion spectroscopy H, i] based on transient processes, 
this approach would have an advantage of being a well- 
quantified steady-state technique, that can be immedi- 
ately coupled to frequency combs and optical frequency 
standards for extreme accuracy and resolution. 

In this paper, we present a numerical and analytical 
treatment of the effect of a narrowband absorption on 
a femtosecond pulse, considered as a dissipative soliton. 
Such a treatment covers both, passively modelocked ul- 
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FIG. 1: Output spectrum of a 100-fs CnZnSe oscillator (black 
solid line) when operated at open air. The atmospheric trans- 
mission (gray) is calculated from HITRAN database [3] and 
corresponds to a full round-trip. The lower graph (b) shows 
the expanded central part of the spectrum. Asterisk denotes 
the absorption line, which is used for quantitative evaluation 
in the last section. 



trashort pulse oscillators with intracavity absorbers, and 
soliton propagation in fibers with impurities. The the- 
oretical results are compared with the experiment for a 
femtosecond Cr:ZnSe oscillator operating at normal at- 
mospheric conditions. We prove that the spectral mod- 
ulation imposed by a narrowband absorption indeed ac- 
curately follows the associated index of refraction when 
the absorber linewidth is sufficiently narrow. 
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2. THE MODEL 

Our approach is based on the treatment of an ultra- 
short pulse as one-dimensional dissipative soliton of the 
nonlinear complex Ginzburg-Landau equation (CGLE) 
d, . This equation has such a wide horizon of appli- 
cation that the concept of "the world of the Ginzburg- 
Landau equation" has become broadly established [lCj . 
In particular, such a model describes a pulse with the 
duration To inside an oscillator or propagating along a 
nonlinear fiber. 

To obey the CGLE, the electromagnetic field with the 



J 



amplitude A(z, t) should satisfy the slowly- varying ampli- 
tude approximation, provided by the relation loq ^> I /To, 
where wo is the field carrier frequency, t is the local time, 
and z is the propagation coordinate. This approximation 
is well satisfied even for pulses of nearly single optical cy- 
cle duration [ll|. When additionally we can neglect the 
field variation along the cavity round-trip or the varia- 
tion of material parameters along a fiber, as well as the 
contribution of higher-order dispersions, the amplitude 
dynamics can be described on the basis of the general- 
ized CGLE ©, ElO 



dA(z,t) 
dz 



[-a + t[P(z, t)] - i 1 P{z, t)}A(z, t)+(a + i^j §^A(z, t)+f [A(z, t)} , 



(1) 



r 



where P = \A\ is the instant field power and a is the 
inverse gain bandwidth squared. The nonlinear terms 
in Eq. ([T]) describe i) saturable self-amplitude modu- 
lation (SAM) with nonlinear gain defined by the non- 
linear operator S, and ii) self-phase modulation (SPM), 
defined by the parameter 7. For a laser oscillator, 7 = 
Aim^lcryst/ )^oA e f f . Here Ao is the wavelength; n and 
712 are the linear and nonlinear refractive indexes of an 
active medium, respectively; l cr yst is the length of the ac- 
tive medium; A e f f — irw 2 is the effective area of a Gaus- 
sian mode with the radius w inside the active medium. 
The propagation coordinate z is naturally normalized to 
the cavity length, i.e. z becomes the cavity round-trip 
number. For a fiber propagation, 7 = 2-Knn2/XoA e ff, 
where n and n% are the linear and nonlinear refractive 
indexes of a fiber, respectively, and A e ff is the effective 
mode area of the fiber [l5[ . Finally, P2 is the round- trip 
net group delay dispersion (GDD) for an oscillator or 
the group velocity dispersion parameter for a fiber with 
P2 < corresponding to anomalous dispersion. 

The typical explicit expressions for S[P] in the 
case, when the SAM response is instantaneous, are 
i) £[-P] = kP (cubic nonlinear gain), ii) = 
k (P — CP 2 ) (cubic-quintic nonlinear gain), and iii) 
p\P] = kP/ (1 + CP) (perfectly saturable nonlinear gain) 
[H LL6(. The second case corresponds to an oscillator 
mode-locked by the Kerr-lensing [8( . The third case rep- 
resents, for instance, a response of a semiconductor sat- 
urable absorber, when To exceeds its excitation relax- 
ation time [T^ . However, if the latter condition is not 
satisfied, one has to add an ordinary differential equation 
for the SAM and Eq. fl} becomes an integro-differential 
equation (see below). 

The (T-term is the saturated net-loss at the carrier 
frequency ujq, which is the reference frequency in the 
model. This term is energy-dependent: the pulse energy 
E(z) = f_ P(z, t)dt can be expanded in the vicinity of 



threshold value a = as a ss 5{E/E* - 1) Qjf], where 
5 — t 2 /go {£ is the frequency-independent loss and go is 
the small-signal gain, both for the round-trip) and E* is 
the round-trip continuous-wave energy equal to the av- 
erage power multiplied by the cavity period. 

The operator T describes an effect of the frequency- 
dependent losses, which can be attributed to an absorp- 
tion within the dissipative soliton spectrum. That can 
be caused, for instance, by the gases filling an oscilla- 
tor cavity or the fiber impurities for a fiber oscillator. 
Within the framework of this study, we neglect the ef- 
fects of loss saturation and let T be linear with respect 
to A(z,t). The expression for T[A(z,t)] is more conve- 
nient to describe in the Fourier domain, A(z,u>) being 
the Fourier image of A(z, t). If the losses result from the 
I independent homogeneously broadened lines centered 
at u>i (relative to ujq) with linewidths Qi and absorption 
coefficients e/ < 0, then the action of operator T can be 
written down in the form of a superposition of causal 
Lorentz profiles [Til H(| 



.4 



i(u - uji)/Qi 



A{z,w). (2) 



In the more general case the causal Voigt profile has to 



be used for T 



A 



[21|. Causality of the complex profile 

of Eq. @ demonstrates itself in the time domain, where 
one has 



1 

f [A{z,t)\ oc^eA J e^ ( ° ! "^ ) ( t - t ')A(z,t')^'- (3) 



The conventional analysis of perturbed soliton prop- 
agation includes approximation of the effective group 
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TABLE I: Laser simulation parameters. The numbers cor- 
respond to a Cr:ZnSe femtosecond oscillator of Fig. [T]with 
Icryst =0.4 cm, w =80 /im, Ao =2.5 /im, n = 2.44, 712 = 10~ 14 
cm 2 /W, ^=0.075, go =2M. 

delay dispersion of the perturbation as Taylor series 
f3'(u>) — (3' 2 {lo - loq) 2 /2 + /3' 3 (lu - lo ) 3 /6, . . . assuming 
that the additional terms f3' 2 , /3' 3 , ■ ■ ■ are sufficiently small. 
This approach is absolutely not applicable in our case, be- 
cause the dispersion, associated with a narrow linewidth 
absorber can be extremely large. For example, an atmo- 
spheric line with a typical width Q = 3GHz and peak 
absorption of only 10~ 3 produces the group delay dis- 
persion modulation of /3' 2 = ±0.9 ns 2 , far exceeding the 
typical intracavity values of (3 2 ~ 10 2 ..10 4 fs 2 . Moreover, 
decreasing the linewidth SI (and thus reducing the overall 
absorption of the line) causes the group delay dispersion 
term to diverge as f}' 2 oc £1~ 2 . 

In the following we shall therefore start with a numeri- 
cal analysis to establish the applicability and stability of 
the model, and then present a novel analysis technique 
based on integral representation in the spectral domain. 
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FIG. 2: (Color online) Part of the pulse spectrum after: (a) 
linear propagation for 25 dispersion lengths inside a fiber with 
two absorption lines and (b) perturbed soliton propagation for 
100 dispersion lengths. Gray curve in (6) corresponds to the 
contribution of only 5R(F) in Eq. ([2]), black curve corresponds 
to the contribution of the complex profile ei = £2 = —0.1, 
Qi = Q 2 = 40 GHz. The initial pulse profile is A(t) = 
Aosech (t/T ), where T =57 fs. 



3. NUMERICAL ANALYSIS 

Without introducing any additional assumptions, we 
have solved the Eqs. (Ill2j) numerically by the sym- 
metrized split-step Fourier method. To provide high 
spectral resolution, the simulation local time window 
contains 2 22 points with the mesh interval 2.5 fs. The 
simulation parameters for the cubic-quintic version of 
Eq. ([1]) are presented in Table |U The GDD param- 
eter P2 = —1600 fs 2 provides a stable single pulse with 
FWHM salOO fs. The single low-power seeding pulse con- 
verges to a steady-state solution during z «5000. 

The pulse propagation within a linear medium (e.g. an 
absorbing gas outside an oscillator, a passive fiber con- 
taining some impurities, or microstructured fiber filled 
with a gas) is described by Eqs. (|1I2[) with zero a, 7, S, 
and the initial A (0, t) corresponding to output oscillator 
pulse. The obvious effect of the absorption lines on a 
pulse spectrum are the dips at uii (Fig. that simply 
follow the Beer's law. This regime allows using the ultra- 
short pulse for conventional absorption spectroscopy 
The nonzero real part of an absorber permittivity (i.e. 

3 ^ ^0) does significantly change the pulse in time 

domain [22j], but does not alter the spectrum. The pulse 
spectrum reveals only imaginary part of an absorber per- 
mittivity (i.e. K(r) ^0). 

Introducing the nonzero SPM coefficient 7 = 
\/32\/AqTq with zero a and S transforms Eq. ([I]) to a 
perturbed nonlinear Schrodinger equation and results in 



the true perturbed soliton propagation. In this case, as 
shown in Fig. [5Jd, the situation becomes dramatically 
different. Besides the dips in the spectrum shown by the 
gray curve corresponding to a contribution of 9?(r) only, 
there is a pronounced contribution from the phase change 
induced by the dispersion of absorption lines (solid curve 
in Fig. [5Jd corresponds to the complex profile of T in 
Eq. ([2])). As a result, the spectral profile has the sharp 
bends with the maximum on the low-frequency side and 
the minimum on the high-frequency side of the corre- 
sponding absorption line. At the same time, the dips in 
the spectrum due to absorption are strongly suppressed. 
In addition to spectral features, the soliton decays, ac- 
quires a slight shift towards the higher frequencies, and 
its spectrum gets narrower due to the energy loss. 

The soliton spectrum reveals in this case the real part 
of the absorber permittivity. However, the continuous 
change of the soliton shape due to energy decay renders 
the problem as a non-steady state case. The situation 
becomes different in a laser oscillator, where pumping 
provides a constant energy flow to compensate the ab- 
sorption loss. 

Let us consider the steady-state intra-cavity narrow- 
band absorption inside a passively modelocked femtosec- 
ond oscillator, where the pulse is controlled by the SPM 
and the SAM, which is described by the cubic-quintic S 
in Eq. Jl]) modeling the Kerr-lens mode-locking mech- 
anism 12]. Such an oscillator can operate both, in 
the negative dispersion regime [l2| with an chirped-free 
soliton-like pulse, and in the positive dispersion regime 
[181 ] . where the propagating pulse acquires strong posi- 



4 




-30 -20 -10 10 20 
co (GHz) 



/(ns) 



FIG. 3: (Color online) Dissipative soliton in an oscillator: (a) 
central part of the spectrum and (b) power P(t). An oscillator 
is filled with an absorber described by Eq. ([2]) with the triplet 
of lines: t\ = 62 = £3 = —0.005, fii = Q2 — = 2 GHz, and 

LJl = -10, LJ2 = 0, LOS = 10 GHz. 



tive chirp. In this study, we consider only the negative 
dispersion regime, the positive dispersion regime will be 
a subject of following studies. 

The results of the simulation are shown in Figs. [3] and 
HI and they demonstrate the same dispersion-like modu- 
lation of the pulse spectrum. Fig. [3] demonstrates action 
by three narrow (17 =2 GHz) absorption lines centered 
at -10, and 10 GHz in the neighborhood of 10 = 0. 
One can see (Fig. that the absorption lines do not 
cause spectral dips at u>i, but produce sharp bends, very 
much like the case of the true perturbed Schrodinger soli- 
ton considered before. One can also clearly see the col- 
lective redistribution of spectral power from higher- to 
lower-frequencies, which enhances local spectral asym- 
metry (Fig. [3^,). Such an asymmetry suggests that the 
dominating contribution to a soliton perturbation results 
from the real part of an absorber permittivity, which, in 
particular, causes the time-asymmetry of perturbation in 
the time domain. This asymmetry is seen in time domain 
as a ns-long modulated exponential precursor in Fig. [3Jd. 

The simulated effect of a single narrow absorption line 
centered at uj = is shown in Fig. [4] for different values 
of peak absorption e and width 17. In Fig. QJi, e = —0.05 
and 17 = 4 GHz (solid curve, open circles and crosses) 
and 1 GHz (dashed curve, open squares and triangles). 
The solid and dashed curves demonstrate the action of 
complex profile @, whereas circles (squares) and crosses 
(triangles) demonstrate the separate action of 3(F) and 
Jft(r), respectively. One can see, that the profile of per- 
turbed spectrum traces that formed by only 3(F) (i.e. 
it traces the real part of an absorber permittivity). One 
can say, that the pure phase effect (3(F), circles and 
squares in Figs. Hk,b) strongly dominates over the pure 
absorption (5R(F), crosses in Figs. Hk,b), like that for the 
Schrodinger soliton. Such a domination enhances with a 
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FIG. 4: (Color online) Central parts of the dissipative soliton 
spectra in an oscillator with the single absorption line centered 
at oj =0. (o) e = —0.05; 17 = 4 GHz (solid curve, open 
circles and crosses) and 1 GHz (dashed curve, open squares 
and triangles). (6) e = -0.005; fl = 4 GHz. The solid and 
dashed curves correspond to the contribution of T described 
by Eq. ((2]). Open circles and squares describe the separate 
contribution of 9 ( T ] . Crosses and triangles describe the 



separate contribution of 5ft ( f 



lowered e (Fig. 0b, crosses), however the |e| growth in- 
creases the relative contribution of 3?(r) and causes the 
frequency downshift of the bend (Fig. 0k). Amplitude of 
the bend traces the e value, while its width is defined by 
17. 

The SAM considered above is modeled by the cubic- 
quintic nonlinear term £ in Eq.([T|). Such a SAM is typ- 
ically realized by using the self-focusing inside an ac- 
tive medium (Kerr-lens modelocking). For reliable self- 
starting operation of mode-locked oscillator it is often de- 
sirable to use a suitable saturable absorber (SA), e.g. a 
semiconductor-based SES AM [2^, [24| • Such an absorber 
can be described in the simplest case by a single-lifetime 
two-level model, giving the time-dependent loss coeffi- 
cient A (t) as 



dA (t) _ A - A (t) 
dt ~ T s 



A(t) 



P(t) 



J S A 



eff 



(4) 



where Ao is the loss coefficient for a small signal, A e ff 
is the effective beam area on the SA, T s and J s are the 
SA relaxation time and the saturation energy fluency, 
respectively. Eq. (|4]) supplements Eq. (Q} and the SAM 
term S in the latter has to be replaced by (Ao — A (£)). 
When the pulse width is longer than the SA relaxation 
time, one can replace Eq. (Q} by its adiabatic solution so 
that 



l + £P(i)' 



(5) 
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where £ = T a J J S S is the inverse saturation power. 

We have simulated Eqs. (|U2l4p in the case of J s =50 fi 
J/cm 2 and T s =0.5 ps, that correspond to the measure- 
ment in Fig. [TJ Two cases have been considered: weak 
focusing (A e ff =4000 fim 2 or saturation energy E s =2 nJ) 
and hard focusing (^4 e // = 1000 /im 2 or saturation energy 
E s =0.5 nJ). We also considered Eq. §5§ for the same 
peak saturation level as weakly focused SA (i.e. £ -1 = 
4 kW). In the latter case the SA effectively becomes in- 
stantaneous, and the perturbed soliton spectrum is the 
same as for the Kerr-lens modelocking, i.e. as for the 
cubic-quintic S. When the saturation energy is suffi- 
ciently large, there is no difference between the models 
expressed by Eqs. ((4]) and (0. The effect of a narrow 
absorption line is similar to that of the soliton of the 
cubic-quintic Eq. (fTJ). Decrease of the saturation energy 
E s causes down-shift of the pulse spectrum as a whole, 
but the narrow bend on the soliton spectrum reproduces 
the real part of the absorber permittivity. One can thus 
conclude that the type of the SAM is irrelevant for an ef- 
fect of the narrowband absorption lines on a dissipative 
soliton spectrum. 

Another important conclusion from the numerical sim- 
ulations is the demonstrated stability of the dissipative 
soliton against perturbations induced by narrowband ab- 
sorption. In the following analytical treatment we shall 
therefore omit the stability analysis. 



The unperturbed soliton parameters are [25 




(G) 



where the equation parameters are confined 



cry 



(7) 



(T > 0. 



Hence, the soliton wavenumber is q = —ja/n. 

Its is reasonable to treat the soliton of the reduced Eq. 
(fTJ) as the Schrodinger soliton with the parameters con- 
strained by the dissipative terms a, a and n (see Eqs. 
(|6I7[) ). This implies that the equation, which has to 
be linearized with respect to a small perturbation co- 
propagating with the soliton without beating, decay or 
growth (i.e. having a wavenumber real and equals to q 
|9|), is the perturbed nonlinear Schrodinger equation 



4. PERTURB ATI VE ANALYSIS OF SOLITON 
SPECTRUM 

To study the transformation of dissipative soliton spec- 
trum under action of narrow absorption lines, we ap- 
ply the perturbation method [§, [25|. Since the ba- 
sic features already become apparent for the perturbed 
Schrodinger soliton and do not depend on SAM details, 
we shall consider the simplest case of the cubic nonlin- 
ear gain = kP. The unperturbed solitonic chirp- 
free solution of such reduced equation with T = is 
a (z, t) = ^osech (t/To) exp [i(f> (t) + iqz] with dtfi/dt = zu. 



^fe^ = 1 YW A{z ' t] ~ %lP{z ' t)A{z ' t) + t [A[z ' t)] ■ 

(8) 

Linearization of the latter with respect to a perturbation 
/ (t) exp {iqz) results in 



2\a(t)\ 2 f(t) + a(t) 2 r (t) 



+T(a + f). 
In the spectral domain, Eq. ((9J becomes 



(9) 



[k (w) - q] f (u) + - [ du' U (u - oj'^J /(w) + — / dJu{u}-J^f*{J^ =S(oj), 



(10) 



where [2f 



U (w) = -T^T^wcsch {ttTquj/2) , 
" IWJ CQsh(7rTo«/2)^ ei l + ( w _ Wl ) 2 /fi a ' [LL) 



6 



Here k (lu) is the frequency-dependent complex wave 
number, and S (lu) is the perturbation source term for 
F corresponding to Eq. ((2]) . 

Further, one may assume the phase matching 
between the soliton and its perturbation. This 
assumption in combination with the equality 
U (lu) = U* (lu), which holds for the Schrddinger 



soliton, results in dco U yu> — lu J f* [lu j = 

The equation (ITU1) for the Fourier image of perturba- 
tion is the Fredholm equation of second kind. Its solution 
can be obtained by the Neumann series method so that 
the iterative solution becomes [23 



fn (w) 



S(m) 3 

k (lu) — q 2ir[k (lu) — q] 



duu U [lu — lu J / n _ 



(12) 



where /„ (a;) is the n-th iteration and fo (lu) — 
S(u)/[k(u)-q]. 

The "phase character" of a soliton perturbation (i- 
multiplier in lhs. of Eq. © and the expression for the 
source term (fTTj) ) demonstrate that the real part of ab- 
sorber permittivity contributes to the real part of soliton 
spectral amplitude. Simultaneously, the resonant condi- 
tion k (lu) — q =0, which is responsible for a dispersive 
wave generation caused by, for instance, the higher-order 
dispersions [9], is not reachable in our case. The reso- 
nance can appear in case of large |e|, k/j, and O/Xb, but 
such regimes are beyond the scope of this work. 

Eq. (fl"2"|) can be solved numerically. Fig. [5] shows 

(fi (u)j (dashed curve) and 3 ^/i (uj)j (dotted curve). 
One can see, that the real part of absorber permittivity 
defines (j (Lu)j while the imaginary part of absorber 

permittivity defines 3 (f (lu)j . That agrees with the 

simulation results and is opposite to the case of a lin- 
ear pulse propagation. One can also see a tiny frequency 
down-shift 9 of the of / (lu) minimum from lui like that 
in the simulations. 

The pulse spectrum (solid curve in Fig. [SJ, results 
from interference of the perturbation with the soliton. 
For the chosen parameters of the absorption line the zero- 
order approximation /o (lu) (open squares) is very close to 
the first-order approximation (solid curve) but is slightly 
down-shifted in the vicinity of the bend maximum and 
minimum. 

With even narrower line width ft of 1 GHz (Fig. [5]) the 
spectral perturbation gets very close to the real part of 
absorber permittivity and, simultaneously, the spectral 

down-shift 9 (location of the 3 f/i) minimum, dashed 

curve) vanishes. The / (lu) (gray solid curves) now per- 
fectly matches the /i (lu) (open circles and crosses) within 
a broad range of e (gray solid curves 1 and 2 as well as 
circles and crosses belong to e =-0.005 and -0.05, respec- 
tively). The bend amplitudes are in agreement with Eq. 

CD- 



A superposition of three identical absorption lines, 
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FIG. 5: (Color online) Central part of the dissipative soli- 
ton spectra perturbed by /i (oj) (solid curve) or /o (lu) (open 
squares) as well as the profiles of 5ft f/i) (dashed curve) and 

(dotted curve) from Eq. (|12|) . Single absorption line 
is centered at lu =0; e = -0.005; = 4 GHz. 

which corresponds to the numerical spectra in Fig. [5J 
is shown in Fig. [7] One can see, that the lowest-order 
analytical solution fo (lu) accurately reproduces the nu- 
merical result. It is important, that a cumulative contri- 
bution of lines into k (lu) does not distort a superposition 
contribution of S (lu) into a soliton spectrum (see Eqs. 
(|llll2l) b This means that the individual contribution of 
a single line within a group is easily distinguishable and 
can be quantitatively assessed, opening way for interest- 
ing spectroscopic applications. 

As Figs. [S] and [5] suggest, the zero-order approxima- 
tion /o (lu) — S (lu) I [k (lu) — q] is quite accurate for a 
description of perturbation in the limit of |e| 1. This 
allows expressing the perturbed spectrum of an isolated 
line (see (|11I12[) ') in analytical form [25| : 
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P(u) = 5(w) + / (w) 



A 2 ^ 2 T 2 sech(^) 2 (^w 2 



1 + n? 



n, 



r 



'/ ) i ir^ 2 



(13) 



Eq. (|T3]) allows further simplification in the case of |e| <C 1 
I 



P(w) « Ao7r 2 sech 




7 (aw 2 + (l + ^9^) ft. 



(14) 



r 



where Eqs. (|6|) and the condition ail 2 <C 1 have been 
used. 

Eq. (|14[) demonstrates that the spectral bend follows 
the real part of absorber permittivity. Spectral down- 
shift of bend is the effect of O (e 2 ) and is not included 
in (|14[) . The perturbation is represented by the term 
in square brackets and its relative amplitude is propor- 
tional to e. Furthermore, the aspect ratio of the kink 
grows with i) the increase of the relative contribution of 
the SAM k/7 ; ii) the gain bandwidth 1/y/a; hi) ap- 
proaching of the resonance frequency w/ to the center of 
soliton spectrum (but the ratio of the aspect ratio to the 
local soliton spectral power increases with \u)A, because 



the former decreases as 



while the latter falls faster 



as cosh (itTouji/2) ); and iv) with an approaching to the 
soliton stability border, which corresponds to vanishing 
a. It should be noted, that smaller a entails the soliton 
width growth (Eq. ©). 

Since the soliton parameters are interrelated, it is in- 
structive to express a through the observable parameters 
such as soliton energy E or the soliton width To. When 
aojf <C cr (e.g. w; «0 or/and an oscillator operates far 
from the stability border a =0), the perturbation ampli- 
tude is inversely proportional to jkE 2 : 



2eiK 



7 (aw 2 



32e ; a 
jkE 2 



7a 



(15) 



For a fixed gain bandwidth, the amplitude scales with 
squared pulsewidth Tq. Ultimately, the latter equation 
is equivalent to 



2e/K 



7 (aw 2 



2f£ 

q 



(16) 



i.e. the relative perturbation amplitude near soliton cen- 
tral frequency is the ratio of incurred loss coefficient to 
the soliton wavenumber, regardless of the z coordinate 



normalization. Therefore, this analytical expression that 
has been derived for the self-consistent oscillator, should 
also be valid for the case of a soliton propagation in a 
long fiber when the conditions of applicability |e| <C 1, 
ill <§C 1/Tq are met. The final form of the soliton spec- 
trum thus becomes 



P(uj) = P (uj) 1 - - > ei 2 

9T 1 + ( "~^ 



(17) 



where Pq(uj) is the spectrum of an unperturbed soliton. 



5. DISCUSSION 

In the analysis above we have shown that, for the case 
of sufficiently sparse, narrow and weak Lorentzian ab- 
sorber lines, their spectral signatures are equivalent to 
the dispersion-like modulation with a relative amplitude 
equal to the peak absorption coefficient over oscillator 
round-trip (or nonlinear length for passive propagation) 
divided by the soliton wavenumber. For quantitative 
comparison with the experiment we recall equation ([5]) 
and express the maximum spectrum deviation of a single 
line at \u — u>i\ = f2; through observable parameters: 



Xi L 



XiL 



0.0319 

IMA*) 2 



(18) 



where xi = 2e; is the peak absorption coefficient of 
the line, L is the absorber path length and Av is the 
full width at half maximum of the soliton spectrum. 
Substituting the actual values of the setup in Fig. Q] 
(/3 2 = -82 ± 40 fs 2 , Av =113 cm-^ 3.39 THz, round- 
trip air path length L = 149 cm, relative humidity 
50 ± 1% at 21 ± 0.5 °C) and taking e.g. the line at 
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co (GHz) 

FIG. 6: (Color online) Central part of the dissipative soli- 
ton spectra perturbed by fi (lj) (open circles and crosses) or 
fo (co) (solid gray curves) as well as the profiles of K (Z 1 ) 

(solid black curve) and 9 f/ 1 ) (dashed curve) from Eq. (|12|l . 
Single absorption line with fi = 1 GHz is centered at co =0; 
e = —0.005 (black solid and dashed curves as well as open 
circles and solid gray curve 1) and e = —0.05 (crosses and 
solid gray curve 2). 




co (GHz) 

FIG. 7: (Color online) Central part of the dissipative soli- 
ton spectra perturbed by /o (co) (black solid curve) and the 
profiles of 5R (dashed curve) and 9 (/°) (dotted curve) 
from Eq. (|12p . Triplet of absorption lines is centered at 
oji = —10, u)2 —Q, and CO3 =10 GHz; ei = £2 = £3 = —0.005, 
fii = fi 2 = fi 3 = 2 GHz. 

4088 cm -1 (122.56 THz, marked with an asterisk), we 
obtain \e/q\ = 0.33 ± 0.02 for the maximum modulation, 
which is in perfect agreement with the observed value of 
31.5% (Fig. [T|o, black line). The agreement is remark- 
ably accurate given the less than optimal resolution of 
the spectrometer (0.25 cm -1 ) and significant third-order 
dispersion of about +10 4 fs 3 that was not accounted for 
in the presented analysis. 



It is important to notice, that the expression (fT8|) 
includes only the externally observable soliton band- 
width and relatively stable dispersion parameter. The 
alignment-sensitive values like saturated losses a, nonlin- 
earity 7, nonlinearity saturation parameter k, etc. which 
are in practice not known with sufficient accuracy, are all 
accounted for by the self-consistent soliton parameters. 

Another important point is the fact, that the signal 
amplitude 2|e/g| can be much bigger than that from con- 
ventional absorption spectroscopy \\L of the cell with the 
same length. The signal enhancement factor can be con- 
trolled by the pulse parameters and it exceeds an order of 
magnitude for the presented case (xiL = 5% for selected 
line and a single-pass cell of a resonator size). For ad- 
ditional sensitivity improvement one can apply the well- 
developed intracavity multi-pass cell technique (26j . The 
expression (IT8|) suggests, that ultimate sensitivity can be 
obtained at the expense of the reduced bandwidth cov- 
erage Av. In this respect, the presented technique has 
the same quadratic dependence of sensitivity on spec- 
tral bandwidth as the conventional intracavity absorp- 
tion spectroscopy Q. 

Further refinement of the presented theory should in- 
clude demonstration of its applicability to arbitrarily 
shaped absorption features. The superposition property 
provides a strong argument for such extension, but it has 
to be rigorously proven for Doppler- and more general 
Voigt-shaped lines, and also for the dense line groups in 
e.g. Q-branches. It would be interesting also to extend 
the theory to the absorber lines at the soliton wings (Fig. 

With the above issues resolved, the soliton-based spec- 
troscopy may become a powerful tool for high-resolution, 
high-sensitivity spectroscopy and sensing. Possible im- 
plementation include soliton propagation in a gas-filled 
holey fibers, as well as already presented intracavity spec- 
troscopy with femtosecond oscillators. The latter, being 
a natural frequency comb source, allows direct locking to 
optical frequency standards, providing for ultimate reso- 
lution and spectral accuracy. 



6. CONCLUSION 

We have been able to derive an analytical solution 
to the problem of an one-dimensional optical dissipative 
soliton propagating in a medium with narrowband ab- 
sorption lines. We predict appearance of spectral mod- 
ulation that follows the associated index of refraction 
rather than absorption profile. The novel perturbation 
analysis technique is based on integral representation in 
the spectral domain and is insensitive to the diverging 
differential terms, inherent to the Taylor series represen- 
tation of the narrow spectral lines. 

The model is applicable to a conventional soliton prop- 
agation and to a passively modelocked laser with intra- 
cavity absorber, the only difference being the character- 
istic propagation distance (dispersion length and cavity 
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round-trip, respectively). In the latter case the prediction 
has been confirmed for a case of water vapour absorption 
lines in a mid-IR Cr:ZnSe oscillator. The model provides 
very good qualitative and quantitative agreement with 
experimental observations, opening a way to metrological 
and spectroscopical applications of the novel technique, 
which can provide a significant (order of magnitude and 
more) enhancement of the signal over conventional ab- 
sorption for the same cell length. 
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